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We study the cosmological stability of a class of theories with a dynamical preferred frame. For a 
range of actions, we find cosmological solutions which are compatible with observations of the recent 
history of the Universe: a matter dominated era followed by accelerated expansion. We then study 
the evolution of linear perturbations on these backgrounds and find conditions on the parameters 
of the theory which allow for the growth of structure sourced by the new degrees of freedom. 
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I. INTRODUCTION 

Theories of modified gravity can be invoked as viable 
solutions to the problem of missing mass. Though ini- 
tially restricted to making predictions regarding qua- 
sistatic, weak field configurations [H, H, Q, more recent 
models are derivable from generally covariant actions 
0, & EJ 0] and so can, in principle, be compared with 
traditional tests of general relativity (for instance see 

The best known proposal, the Tensor- Vector-Scalar 
(TeVeS) theory of gravity pro posed by Bekenstein has 
been studied in some depth [ill . [l2l [T3L flij , and has been 
found to explain a plethora of observations, from galactic 
rotation curves to the growth of structure during recom- 
bination without invoking dark matter. We have shown 
[IH that it can be seen as a particular case of Einstein- 
Aether theories of gravity [l6| with non-canonical kinetic 
terms and that one can consider a different subclass of 
these models which seem to be observationally viable [I?) • 
For convenience and to distinguish them from standard 
Einstein- Aether theories, we have named these theories 
Generalized Einstein Aether Theories. 

In this paper we wish to explore, in as general terms 
as possible, the cosmological stability of these theories. 
It has been shown that the "vector" part of TeVeS plays 
an essential role in the growth of structure that leads to 
the formation of galaxies [l9[ . Perturbations in the time- 
like vector field are gravitationally unstable and grow in 
a manner akin, but not identical, to that of normal, pres- 
sureless matter. We wish to find if this phenomenon is 
present in Generalized Einstein Aether Theories. 

In this paper we present the following results: 

• we present the full set of background, and linear 
perturbation, equations for Generalized Einstein 
Aether Theories; 

• for a broad class of models there is accelerated ex- 
pansion following a period of matter domination; 

• in a pure baryonic universe, the growth of structure 
can be sustained and enhanced by the presence of 
the Aether vector field; 

• although the growth of structure is driven by the 
vector field, unlike in TeVeS, it remains strictly unit 



timclike leading to qualitatively different features; 

• there are choices of parameters in the Einstein- 
Aether theories that give a reasonable fit to current 
measurements of large scale structure 

• a few notable features can be used to distinguish 
it from conventional gravity and dark matter, in 
particular a different growth rate for structure and 
the presence of non zero, scalar anisotropic stresses 

There are a number of steps we must go through be- 
fore we can assess if there is the desired instability of the 
vector field, or aether, in these theories. We are look- 
ing for it in the cosmological setting and hence we must 
identify the restricted class of theories that give us ac- 
ceptable cosmologies. We do so in Section IPlTl Once we 
have established what region of parameter space we want 
to explore, we then check if perturbations in the vector 
field grow in a matter dominated era. We shall see, in 
Section [V] that this defines a subset of parameters. We 
then consider the full coupled system of metric, matter, 
and vector field perturbations for modes outside the hori- 
zon well before recombination and trace their evolution 
up to the present time. 

As will become clear, we are able to identify the subset 
of theories which could possibly lead to realistic growth of 
structure in a purely baryonic universe. In Section IXl we 
discuss if the constraints we find in this paper are compat- 
ible with a host of constraints that have been identified 
in [13 . 



II. THE THEORY 

Let us first recap the structure of Generalized Einstein 
Aether Theories. A general action for a vector field, 
coupled to gravity can be written in the form 



S = 



hl 2 PL R + C(g^,An 



S 



M 



(1) 



where g^, is the metric, R is the Ricci scalar of that met- 
ric, Mpl = , 1 is the Planck mass, Sm is the matter 



action, and C is constructed to be generally covariant 
and local. Sm couples only to the metric, g^ v and not 



2 



to . Throughout this paper, repeated Greek indices 
are understood to be summed over 0..3 and Latin indices 
are summed over 1..3. Furthermore we will consider for 
most of this paper that is unit time-like. 

In this paper we will focus on the subset of theories 
proposed in [171 ] . i.e. we consider a Lagrangian that only 
depends on covariant derivatives of A^ 1 and the time-like 
constraint. It can be written in the form 

£(g ab ,A a ) = ±M 2 M 2 PL T(IC) + hl 2 PL \(A a A a + l), (2) 



where 



K = 

-a/3 _ 



M- 2 K, a \ a V a AiV p A°, 



Here Ci are dimensionless constants and M has the 
dimension of mass. A is a non-dynamical Lagrange- 
multiplier field with dimensions of mass-squared. One 
may additionally consider a l C4 term in K a p 1IJ propor- 
tional to ga/jA^A 1 ' [Tq . In [ijj we considered the limit 
where matter sources in the gravitational field equations 
are negligible. We found that for the perturbed aether 
degrees of freedom to have a positive Hamiltonian, and 
to avoid superluminal propagation in spin-0 aether per- 
turbations and gravitational waves collectively implied 
that: 



ci < 0, c 2 < 0, and c x + c 2 + c 3 < (3) 
The gravitational field equations are 

G a() = T a/3 + 8irGT™p atter (4) 

where the stress-energy tensor for the vector field is given 
by 

1. 



where 



dT 



dK. 1 



and 



Y af ) = ci[V s A a V(:A - V a A 6 V A^ 



Brackets around indices denote symmetrization. 
The equations of motion for the vector field are 



(5) 

(6) 
(7) 

(8) 
(9) 



Variations of A will fix A^A^ = — 1. Transvecting the 
equation of motion for A^ with A 11 allows us to solve for 
the Lagrange multiplier: 



(10) 



III. BACKGROUND EVOLUTION 

We first consider the case of a homogeneous and 
isotropic universe in which the metric is of the form 



g^ydx^dx" 



-dr 



a{t) 2 5ijdx L dx\ 



(11) 



where t is physical time and a(t) is the scale factor. 
These solutions will be used as the background on which 
we can study the growth of perturbations. The vector 
field must respect the spacial homogeneity and isotropy 
of the system and so will only have a non-vanishing l f 
component; the constraint fixes A^ — (1,0,0,0). The 
energy-momentum tensor of the matter is of the form 



T -matter 
a/3 



pU a U p + P(g a p + U a U p ), 



(12) 



where p is the matter energy density, P is pressure and 
we have introduce a unit timelike four- vector U u ', i.e. 
satisfying gjlv U»U» = -1. 

The equations of motion simplify dramatically with 
these symmetries and we find that: 



K = 3 



aH 2 
M2~ 



(13) 
(14) 



where H = -, the dot denotes differentiation with re- 
spect to t, and, following [Ta |. we define a = c\ +3c 2 + C3. 
Recall that a is negative and hence so is K,. 

The modified Einstein's equations now become: 



[1 - T K a\H 2 
2aT K \R 2 - 2[1 



-TM 2 
6 



P (15) 



2 a 
1 



T K aH - -TM 2 = 8ttGP. (16) 



We can rewrite these equations in a different form: 



H 2 - 



8ttG 



(17) 



dt 



{-2H + T K aH) = 8irG(p + P). (18) 



As discussed in [17| , the form of T for negative values 
of JC cannot be constrained by looking at the structure of 
the modified force law near spherically symmetric bodies 
or near the Newtonian regime where K, > 0. Hence we 
can consider any general form. Initially it is informative 
to restrict ourselves to 



T = l(rK.) n . 



(19) 



With such a choice, we find that equation [T7] is more 
clearly a modified Friedmann equation: 



1 



2(n— 1)' 



H 2 



8ttG 



(20) 
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where 



e = (1 - 2n)7(-3a) n /6. 



(21) 



For the purpose of this paper and as can be seen from 
the Friedmann equation, we have restricted ourselves to 
flat universes. This implies that there is a relationship 
between 



Q m = 8ngp Q /3H 2 (22) 
(where Hq is the Hubble constant today) and 7. Indeed, 

6(n m - 1) 



7 



M 

5o 



2(1-1) 



(1 - 2n)(-3a) 
We can immediately identify a few special cases: 



(23) 



• n = 1/2 - the Friedmann equations are unchanged 
(e = 0) and there is no effect on the background 
cosmology; 

• n = 1 — we have that e = 7a/2 and recover the 
results of fl8] (7 < 0), i.e. Newton's constant is 
rescaled, G' =G/(l + e); 

• n = - we recover a cosmological constant, A ~ 
sign(— ~f)M 2 . 

More generally we will obtain different regimes depend- 
ing on the relative size of each term in equation [^DJ Con- 
sider e < 0. We can define e = — e 2 to rewrite the modi- 
fied Friedmann equation as 













* +i (£) v 









8ttG 



p. (24) 



Clearly H -> i? eg where ff e? = Me 1 /' 1 "™'. We can ex- 
pand around the equilibrium solution, H ~ H eq (l + /i) 
to find the rate of approach to equilibrium given by 



1 



87rGp 



2(1 - n) 3H 2 q 



(25) 



In this regime accelerated expansion is approached as 
p — > 0. Note that we have assumed that e < so that 
we have n < \ and 7 < or n > | and 7 > 0. Note also 
that this corresponds to O < 1. 

There is an important qualitative difference in the evo- 
lution of H depending on the value of n. If n < 1, then 
the exponent of H/M in equation (f2T)]) (and in equation 
(|24[)) is smaller than the exponent of the standard iJ 
term. The H/M terms will therefore contribute most 
when \H/M\ is small. With H currently positive and 
smaller than in the past, this means that we could be 
approaching a phase of constant H. On the other hand, 
if n > 1, then the H/M terms will contribute most when 
\H/M\ is large, in other words in the past. Since the 
expansion is observed to be approaching a period of ap- 
proximately constant H (w ~ — 1), ife < we prefer 
n < 1. 



We can now consider e > 0. For this regime there is 
again a competition between the conventional H 2 term 
and the modification due to the vector field. It is still the 
case that if n > 1 any modifications in the background 
cosmology come into play for high values of H/M, i.e. at 
early times. As claimed in [l7j . this can be a source of 
cosmological inflation in the early universe. This is not 
the regime we are interested in, in this paper. 

For < n < 1 we find that the expansion rate at late 
times will be retarded. If we assume p oc a~ 3 we find that 
a oc t 2 ™/ 3 with the onset of the regime occurring when 
H ~ e -1 / 2 ( n-1 )M. Note that as e -> 0, this will occur 
as t — > 00. Holding 7 and a fixed, this means that as 
n — > 1/2, t — > 00. Also, since H ~ 2n/3 in this regime, 
as n tends to the expansion rate decreases. 

For n = 0we find that the equations correspond to an 
Anti-deSitter spacetime with additional matter; the Uni- 
verse will expand, slow down and then loiter at H ~ 
and constant a. This point triggers the onset of an exotic 
range of behaviour, for n < 0. For negative values of n 
we find that a will expand and H will tend to f^l^M 
after which H changes sign discontinuously and a con- 
tracts. This is signalled by H — > 00 at the transition. 
The energy momentum tensor transits smoothly at this 
point, guaranteeing that the system is well behaved. 

We can summarize the behaviour in a plot in (n, 7) in 
Figure [TJ The region of interest to us is clearly 7 < 
and n < 1/2 and 7 > and n > 1/2. In this region we 
can envisage, for example, a period of matter domination 
after recombination followed by late time accelerated ex- 
pansion. 



IV. LINEAR PERTURBATION THEORY 

We now look at linear perturbations on the homoge- 
neous background. Throughout this paper we will focus 
on scalar perturbations and work in the conformal New- 
tonian gauge. The metric can be expanded as 



(1 + 2e^)dt 2 + a 2 (l - 2e$)6 ij dx i dx j , 



(26) 



while the four vector A a can be expanded as 

A" = (l + eX,e&>Z) (27) 
= (l + eX^djZ). 

Here indices have been lowered with the background met- 
ric. In Fourier space: 



A 1 " = (1 - eV,i-kjV), 



(28) 



where, for computational convenience, we have defined 
V = Z/a and have used the fact that the constraint fixes 

X = 

We can expand K. to linear order: 



(29) 
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FIG. 1: A schematic representation of the types of the late- 
time background evolution as a function of (n, 7) for n < 1. 



K° = 3 



aH 2 
M 2 

M 2 V a 



3i7* + 3$) 



(30) 
(31) 



We now switch to the conformal Newtonian gauge. 
Derivatives with respect to conformal time, 7/, are de- 
noted by ' and — = 7i. The evolution equation for the 
perturbation in the vector field becomes 







d[V" + k 2 V + 2HV' + 2H 2 V 
+ $' + 2Hf] 

a" 



(32) 



c2[k 2 V + QH 2 V - 3—V + 3$' 
a 
11 

2t/ 1 o«2 



ft 



+c3[k 2 V + 2H 2 V V 

a 



Fkk, 



K e aH 



F K 

-K°'(- Cl (V' 



*) + Z02HV + c 3 HV)}. 



The perturbation in the vector field is sourced by the 
two gravitational potentials $ and ^ . These in turn are 
seeded by the vector field and the perturbations in the 
matter field. The first order perturbations to the vector 
field's stress energy tensor are: 



a 2 5f° 



F KCl [-Hk 2 V - k 2 V' - fc 2 *] 
+F K a[2Hk 2 V + 6H$' + 6H 2 <F] 



(33) 



1 



-afF K M 2 K 6 - iF KK aH 2 K' i 



a 2 5f° 



a 2 5f\ 



F K ci[-Uk 2 V - k 2 V' - k 2 ^} 
+F K a{2n - l)[Hk 2 V + + 3H 2 *], 

, : = ihF KCl [V" + 2HV' + —V (34) 

+ nv] 

a" 

+ihF K a[2H 2 V - — V] 

+ik l F KK K Q '[c 1 {UV +Y + $) 
-aHV], 



F K c 2 k 2 [2nV + V']8 i j (35) 
F K a[k 2 HV + + <P" 



+2^ 
a 



2H 2 ^ 



+F K ( Cl + c 3 )[2HV + V']kikj 

+ la 2 F K M 2 K e S\ 
2 J 

+F KK [-aK e — - (cj + c 2 + c^K'H 2 
a 

-uHK e ' + aK°'<P' + 2aK°'ny 
-a\n{F KK )'K e H + c 2 K°'k 2 V]5 l 1 

+(ci + C3)F KK K°'Vk i k j . 

where the second expression for a 2 ST° assumes the 
monomial form for F(IC). In the absence of anisotropic 
stresses in the matter fields, we may obtain an algebraic 
relation between the metric potentials <& and '5 by com- 
puting the transverse, traceless part of the perturbed 
Einstein equations, i.e.: 



k 2 {^ - $) 



- (X k-i k j 



(36) 



= (ci + c 3 )k 2 [F K (2HV + V) + F KK K°'V]. 



We can see then that generally T *j 
From the perturbed Einstein tensor we obtain 



AS'j if ci 



a 2 [G% 



Y2 Hik 3 G° 3 ] 



2k 2 <f>. 



c 3 = 0. 



(37) 



Similarly, 



a 2 [f° - —Hikjf^] = -F KCl k 2 [V + * + (3 + 2c 3 )HVl?8) 

where we have used the identity provided by (|3"2f . 

We take the dominant remaining contributions to be 
from the baryons and radiation, treated as perfect flu- 
ids. No contribution from anisotropic shear is considered. 
Therefore, the perturbations to a component described 
by stress energy tensor T are given by: 



ST° t 



5T' 
ST 



0(a) 
) 

i(a) 



-p(a)S(a) 

\P(a) + P(a)]v l (a) 

[P(a) + 5P{a)W 3 



-8T\ 



(39) 
(40) 
(41) 
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where p( a ) and Pr a \ are the background energy density 
and pressure of the fluid, S/ a ) is the density contrast 



P(a) ' 



SP( a ) is the perturbation to the fluid pressure, and v 1 ^ 
are components of the fluid's 3 velocity. We also define 
the divergence of the fluid velocity, 9 — ikiV 1 . Calculating 
the contribution from this tensor in Einstein's equations, 
we find the following expression for Poisson's equation: 



k 2 <£ 



V. 



--F KCl k 2 [V' + V + (3 + 2c 3 )HV] (42) 



EVOLUTION OF VECTOR 
PERTURBATIONS 



We now investigate solutions to the vector perturba- 
tion equation (|32p . We first establish the conditions for 
the existence of growing modes during the matter domi- 
nated epoch, since in their absence the low amplitude of 
perturbations in the matter at last scattering, and the 
subsequent effects of Silk damping, will result in a failure 
to evolve structure in the universe. 

We consider therefore that we are following the evo- 
lution of the vector field modes during an epoch when 
the scale factor a(rf) oc 77™. to = 2 during the matter- 
dominated epoch. We also suppose that F(fC) is domi- 
nated by a single power law - as above F(JC) = 7 (—fC) n 
- and that this form is consistent with a period of radia- 
tion domination followed by matter domination. In this 
case the source-free part of equation (f32|) is: 



where 







bi 
hi 



V" + b 1 -V + {b 2 + b z {kj 1 f)\v i (43) 



= 2m- (2m + 2) (n- 1) (44) 

= (to 2 + to + 2m(n — 1)(to + f ))q (45) 

+(m 2 — to — 2m(n — l)(m + 1)) 
2 

63 = l + c 2 +c 3 + -(n-l)a. (46) 

Here ci = Cj/ci, and a = a/ ' c\. Because of the condition 
©, c~2 > 0, and f + c 2 + £3 > 0. Consequently a > 
also. 

In addition to the source-free terms included in P3"|) . 
the vector evolution equation (f32|) contains both self- 
gravity terms and terms sourced by matter perturba- 
tions. Since we are interested here in the initial growth of 
structure, early in the matter dominated epoch, we may 
ignore the matter perturbations. What then of the self- 
gravity terms in (f3"2")) ? To lowest order (in e) these may 
be collected into the following instructive form (and di- 
vided by Ci) in which they can be compared to the terms 
of equation (|4*3")l : 



(*' - + (1 + a (2n - I)) ($' + 



(47) 



-2(n- 1) 2a 



1 I 



TO Cl 



(48) 



From ([36]) we see that the first term (vff' — $') is pro- 
portional to Tk = dT I dK and its derivatives with re- 
sepct to K and r\. Recall that T = j(—JC) n (equa- 
tion HI])), and that we have found that interesting back- 
ground cosmology suggests that (at least at late times) 
we confine our attention to n < 1. Equation (|13|) . 



K 



aH 2 
' M 2 ' 



(49) 



and the fact that in order to recover galaxy rotation 
curves we must have M be of order the current Hubble 
parameter H ([ljj]), together imply that unless \a\ <C 1, 
JC will be large at large redshift when H ^> H n . Thus, in 
the period of interest soon after matter domination, we 
expect Fk and its derivatives to be very small, and so 
the first term in (|47|) to be negligible. 

The second term in (j47|) is proportional to ($' + 7i<I/). 
According to the G°j Einstein equation 



k 2 ($' + WS) = -H 2 (l + 3w)0 



(50) 



We are specifically interested in the growth of vector field 
modes at times when matter seeds are small, so we can 
set 9 = 0. Meanwhile ST°j is, according to equation (|39[) . 
proportional to Tk and its derivatives. As we have just 
seen, such terms are suppressed at early times. 

Finally, the third term in equation (|47|) is proportional 
to \P . Using the Poisson equation (|4"2"| , and equation (|36|) , 
one can reexpress ty: 



* 1 



1 



= l -T K {c 1 + 2c z )V' 
+T K U[^ (ci + 2c 3 ) 



(51) 



777 -L- 1 

-2( Cl + c 2 )^—(n-l)]V 



4ttGci 2 p 
P 



[6 + 3 



W 

k 2 J 



Again, all terms on the right hand side are either sup- 
pressed by factors proportional to Tk-, or are sourced by 
matter. 

We thus conclude that for an analysis of the existence 
of vector growing modes soon after the onset of matter 
domination in the absence of matter seeds, we can confine 
ourselves to the source- free equations (|43|) . 

We can solve analytically the source free vector per- 
turbation equation (|43[) even with the inclusion of the 
inhomogeneous terms so long as we confine ourselves to 
a particular cosmological epoch, i.e to fixed to. The so- 
lution is: 



Ml-b! 



-/ 2 (k)F(/3, Vhkrj)}, 



(52) 
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where J and Y are Bessel J and Y functions respectively 
and the fa are functions to be fixed by boundary condi- 
tions. 

We look first for the homogeneous (krj = 0) solution 
to (|4"3")) . If p > 0, the appropriate form of a solution is 
V oc rj p , hence 



(53) 



The Ci are constants of integration and for convenience 
we have defined (3 = \ {1 - 46 2 - 26i + 6f)s. 
A growing mode therefore exists if and only if 



during the radiation era. 

We find decaying solutions of the type (|53"j) for < a < 
1/24 during the matter era and < a < 1/8 during the 
radiation era and that, signalling the transition [3 to an 
imaginary value , greater respective values of a lead to 
decaying oscillatory modes. 

For more general F(JC), only in the regime (n < 1/2) 
do we both find a growing mode, and an approach to a 
late time acceleration in the background. Therefore we 
will restrict our attention to these values of n for the 
remainder of the paper. 



i(l-6i)+/?>0. 



(54) 



VI. EVOLUTION OF COUPLED 
PERTURBATIONS DURING THE MATTER ERA 



In the matter era (m=2), this condition describes three 
regimes where growing modes occur: 



< a < 



a > 



1) 



6n — 7 
3(2n 

6n — 7 
3(2ra- 1) 



for 7/6 < n 
for n<l/2. 



(55) 

(56) 
(57) 



where we have insisted that a > 0. 

In the radiation era (m=l), alternative regimes lead to 
growing modes: 



n - 2n ~ 2 c 
< a < for 1 < n 



a > 



(2n - 1) 
2n - 2 
(2n - 1) 



for n < 1/2. 



(58) 

(59) 
(60) 



If n is unchanged from the radiation-dominated era 
to the matter dominated era, then the requirement of 
growing modes in the matter dominated automatically 
results in growing modes in the radiation-dominated era. 

We first consider the scaling value (n = 1) previously 
well-explored by Jacobson, Mattingly and collaborators. 
Indeed in this case we have that simply that: 



bi = 4 
b 2 = 6a 

1 



P = (4-6S 



(61) 
(62) 

(63) 



during the matter era and 



bi 
b 2 

= 



2 

2q 
1 



(64) 
(65) 

(66) 



We now consider the effect of the vector field during 
matter domination; in doing so we take the dominant 
remaining contribution to the energy density to be bary- 
onic, treated as a pressureless perfect fluid with energy- 
momentum tensor T. 

V a <5T afc = produces the following: 

5' = -6 + m' (67) 
9' = -W + kH. (68) 

Furthermore we introduce the variable: 



V = E 



(69) 



Examining (|32p we expect a first order evolution equa- 
tion for E. Therefore, the perturbation variables are: 
S, 9, VP, V, E. In addition we have two constraints be- 
tween the variables, given by: 

L ^ K cik 2 [E + * + (3 + 25 3 )W](70) 



-F, 



/fc 2 (*-$) = 



-4irGa 2 p[6 + 3H—} 

( Cl +c 3 )k 2 [F K (2HV 
+F KK K°'V]. 



E) 



(71) 



Following [20( we choose to eliminate 4" via the con- 
straint equations and integrate the remaining variables. 
We use the G° component of Einstein's equations as a 
first order evolution equation for $ and reserve Poisson's 
equation as a consistency check. 

We can now examine the evolution of perturbations 
in the matter era after recombination. In solving the 
full system of equations numerically, we have found that 
the solution (|53|) describes the evolution of V during the 
matter era for universes that produce a realistic descrip- 
tion of the matter power spectrum today, vindicating our 
choice to neglect the effective metric and matter sources 
to the vector field equation during the matter era. 

For ease of illustration we will initially consider only 
the case where V is described by a growing monomial, 
i.e. 



V = Vn I - 

Vo 



(72) 
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Taking the vector field to be described by (|72|) , the ar- 
guments from Section [V] may be similarly used to show 
that, concomitantly the components of vector field stress 
energy tensor 5T a b will be dominated by those propor- 
tional to V and E. It follows then that during the matter 
era both the quantity a 2 5T° and the difference between 
the conformal Newtonian gauge potentials take particu- 
larly simple forms: 

-l E £(k)k 2 V 5+p - 6n 

-z s £(fc)fcV +p ~ 6n , 



a 2 ST° n 



o 



(73) 
(74) 
(75) 



where 

m 

Ie 
Is 



lV {k) 



j,6-6n 
'"hub 



-(c 1 (2+p)n + 2a(l 
-(ci + c 3 )n(6n — p - 



3aQ T 

- 2n)n) 
10). 




(77) 
(78) 



Vtoday 



where k hub = 

The vector field affects our evolution equations for the 
matter and metric perturbations only through the above 
two expressions. The first may be interpreted as propor- 
tional to an energy density contribution of the vector field 
(hence Ie) in the G° Einstein equation, proportional to 
both the co-moving wavenumber squared and the confor- 
mal time rj raised to the power 5 + p — 6n. The second 
contribution may be seen as an anisotropic stress due to 
the vector field (hence Is) and has an identical k and r\ 
dependence as the energy density. Indeed, for general c, 
it is clear that the U and thus both terms should be of 
similar order. 

We will now characterize the possible evolution of the 
perturbations during the matter era for a wavenumber k 
that is outside the horizon at the beginning of integration 
{krj « 1) and inside by the end {krj >> 1). Recall 
that in the absence of an vector field, and for suitable 
initial conditions, the potential and density contrast are 
constant outside the horizon, whereas inside the horizon 
the potential remains constant and the density contrast 
grows as rj 2 . In the presence of the vector field then the 
full system of evolution equations in the matter era is as 
follows: 



5' = 



3W$' + fc 2 $ + 3H 2 $ = 



-6> + 3$' 
-H6 + k 2 <P 
+lst{k)k 2 V 5+ P- 6n 
~ATrGpa 2 6 



i 



(79) 
(80) 

(81) 

2 5+p-6n 



terms which necessarily depend on positive powers of k 
with respect to other terms. Therefore the velocity diver- 
gence 9 becomes decoupled from the evolution equations, 
which may then be combined to obtain a second order 
evolution equation for 5: 



5" + -5' = 6(5 + p - 6n)? s e(fc)^ 3+p " 6 " (82) 
V 

with solution: 



S = d(k) + 



6/s£(fc) 



(10 + p-6n) 



5+p— 6n 



(83) 



where C\ is a constant of integration and we have omitted 
the decaying mode. 

Therefore even before horizon crossing, the anisotropic 
stress term due to the vector field can influence time evo- 
lution of the baryon density contrast to a degree dictated 
by the parameters of the theory and initial conditions on 
the vector field. 



B. krj » 1 

In the limit kr\ >> 1 the fc 2 $ term in the quan- 
tity G° dominates the terms with time derivatives, e.g. 
7^$' ^ $ /rj 2 . The system of equations is given to a good 
approximation by: 



6' 



k z § 



-6> + 3$' 
2 



(84) 



6 + k 2 $ - lst(k){kr ] ) 2 ri 3+p - 6n (85) 
V 

^Ettk)(kv)v 3+p - 6n - \& 
2 rj z 



(86) 



Differentiating the density equation and substituting 
in for 6" using the velocity equation, we recover: 



5" = -9- k 2 <P + l s (kri) 2 ri 3+p - 6n + 3$" (87) 
V 

= -1 S ' + -$' - fc 2 $ + Z s (fc?7) V +p - 6n + 3$" 

As before we may neglect the contribution due to time 
derivatives of $, then using the Einstein equation to sub- 
stitute in for fc 2 $ we finally recover: 



5" + -6> -^8= (h E + ls)d(k)(kv) 2 ri 3+p - 6n 
with solution 



A. k-q « 1 

We assume adiabatic initial conditions for the fields 
<5, $, 8. In the super-horizon limit, one we may drop all 



S(k, V ) = C 2 {k)ri 2 



(¥e+Is) 



(5+p-6n)(10 + p-6n) 



(89) 
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FIG. 2: Influence of the vector field on the growth of 8. Dotted 
line is evolution of <5 for a GR/baryon universe. 8 is constant 
until kr\ 1 and then grows as rf . We now include the vector 
field. We choose n = 0.4 and p = 0.4. Dot-dash line shows 
the vector field initially subdominant to the matter pertur- 
bations. For kr\ « 1 it has negligible influence but at later 
times its energy density dominates and leads 8 to grow as 
^7+ P -6n _ ^5 gQijj i me snows a vector field dominant 

even at early times. Before horizon crossing the anisotropic 
stress dominates and leads to 8 growing as jf +p ~ (in — rf . 
Eventually again the energy density dominates and 8 grows 
as -rf (parallel with dot-dash line) 



where C2(k) is a constant of integration. 

We can see then that for sub-horizon modes, the influ- 
ence of the vector field on the evolution of 8 is a combi- 
nation of the effect of the energy density and anisotropic 
stress contributions though both, in this limit, result in 
the same contributions to the scale dependence and time 
evolution of the density contrast. Two examples of the 
influence of the vector field on the evolution of a baryon 
overdensity mode are illustrated in figure (|2"|). 



VII. THE POWER SPECTRUM OF MATTER 
PERTURBATIONS 

Thus far we have explored the dynamics of the various 
instabilities of both the vector field and baryon perturba- 
tions in generalized Einstein- Aether theories. Under cer- 
tain conditions, structure does grow. We can now make a 
first attempt at calculating the power spectrum of baryon 
fluctuations in a realistic universe. We will forego the 
complete analysis using a full Einstcin-Boltzmann code 
and focus on the lower order moments of the matter field. 
Furthermore, we will not include Silk damping during re- 
combination at this stage. Like dark matter, the vector 



field is coupled to the other fields exclusively via source 
terms from the metric potentials and so as in the case 
of dark matter we expect the above approximation to 
be accurate when looking to evaluate the linear baryonic 
matter power spectrum [201 ] . 

The matter density can be split up into a radiation and 
baryonic component, p = ps + p 7 with ps oc a~ 3 and 
p 1 oc a -4 . We set equality pb = p-y at a eq = f2^/f2 m . 
We now split the evolution of perturbation into two 
stages, before and after recombination at a* = 10 -3 . Be- 
fore recombination, the baryons and radiation are tightly 
coupled through Thompson scattering and they can be 
described by a combined fluid: 



S'r 



B~, 
I 

B-, 



9 B7 + 3$' 

n „ 



l + R 



'B-y 



R k 2 

i + i?T 



(90) 

6 B -y + k 2 t>, (91) 



where R = 4p 7 /3pB- 

After recombination, the two fluids decoupled and 
baryons will evolve as in section [V] while the radiation 
will obey: 



= -3^+4$' 



k 2 

9' = 

7 7 4 



'—S^ + k 2 ^. 



At recombination, the matching conditions are: 



3£ 7 /4 



and 



'B — c 7 



'Bf 



(92) 
(93) 

(94) 
(95) 



Therefore our approach is as follows: 

• For a given set of flb^r,H 07 a 7 n, and 7 determine 
the background evolution of the universe. We fix 
h = 0.73 where Hq = 100h(km/ s)Mpc~ 1 as well as 



Q h h 2 = 0.023 and Q r h 2 = 2.47 • 10" 



Choose initial 



conditions 
and E. 



for 
We 



the 
choose 



fields 
adia- 



batic initial conditions for fields other other than 
the vector field which, due to the smallness of !Fk, 
will be decoupled from the other fields at early 
times. 

• During tight coupling, numerically evolve the fields 
and V using the G° Einstein equation and 
the vector field equation of motion while eliminat- 
ing 'J' via the constraint equation (|36p . Then evolve 
the matter fields using the tight coupling condition 
and equations (|9"0|) and Q9ip. 
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At recombination use the matching conditions (|94[) 
and (|9"5|) . Then continue to evolve the vector field 
and metric potentials as before but now evolve <$ 7 
and # 7 using (f9"2")l and (|9"5)) while evolving <5& and Ob 
using ([67)) and (f68|) . We cease integration at z = 0. 



The observable today may be read from (|89|) . It 
contains two functions which are determined in the nu- 
merical evolution: C2(k) and The form of the func- 
tion C 2 (fc) will follow from a given set of parameters and 
initial data. Meanwhile, the function £(fc) is determined 
by the manner in which V reaches a dominant growing 
mode solution. We may isolate limiting, analytic cases 
for the modes of interest: The case in which the initial 
values of E and V/rj are significantly smaller than the ini- 
tial metric source terms (Case 1) and the case in which 
the initial values of E and V/rj are significantly greater 
than the initial metric source terms (Case 2). 



A. Case 1 

We first consider modes that enter the horizon after 
radiation contributions may be neglected. The scale de- 
pendence of these modes may obtained to a good ap- 
proximation by matching our super-horizon solution to 
our sub-horizon solution at horizon crossing. A given 
mode will enter the horizon when krj en t ~ 1 where rj en t 
is the conformal time at horizon crossing. Matching the 
solutions from the previous section we find the following 
solution for C 2 (k): 



C 2 {k) 



fiprim (J$)k 



61, 



■(10 +p- 6n) 



(96) 



(hlE + ls) 



(5 + p — 6n)(10 + p — 6n) 



m)k 



-3— p+6ro 



where S pr i m (k) is the primordial amplitude of S. If we 

assume adiabatic initial conditions and a scale invariant 

-3/2 



spectrum for $ then 8 pr i m oc k~ A l 1 . We must now find 
the scale dependence of To do this it is necessary 

to examine how V reached the monomial solution. To do 
this we solve (p?2"j) for super-horizon modes during radia- 
tion domination, which reduces to: 



V" 



bi 



— V — — * 

9 v ^ prim 

Tj 7j z V 



V 



(97) 



where I3 = — [a(2n — 1) + 5 — 4n] . 

We have assumed that at such early times the vector 
field does not affect superhorizon metric modes enough 
for them to deviate appreciably from their primordial val- 
ues in determining the early evolution of V. The equation 
may be solved to yield: 



V = Dnf+ +D 2 r, p - 



h+b 2 



We now impose V(k)/rj << $> pr im- This implies that 
the decaying mode coefficient D 2 should be very close to 
zero, and it shall be neglected from now on. Now differ- 
entiating the solution, the condition V'(k) = E(k) « 



prim 



requires that 



prim 



{k)vl' 



p+{h + b 2 ) 



(99) 



and so 



V = * 



prim 



(fc)ry 



h + b 2 



(100) 



Therefore, providing that p + — 1 > 0, the complemen- 
tary function solution will tend to overcome the partic- 
ular integral and the monomial solution will be reached. 
Then comparing with [72] we find for these modes that: 



Vo(fc) = ^> pnm {k) 



-m 



(101) 



Thus, from equation (|76|) we see that both Vo(k) 
and £(fc) will have the same k-dependence as ^ pr imik). 
Putting everything together, we find that for modes en- 
tering the horizon after radiation domination that: 



8 2 {k^) = f^k 



(k « k eq ) 



(102) 



where the function /1 (77) is determined by the solution to 
the field equations. 

Therefore for modes that enter the horizon after ra- 
diation domination, the scale dependence of the power 
spectrum is expected to be identical to that of a pure 
baryonic or ACDM universe up to a multiplicative func- 
tion of the conformal time. 

We now look at the behaviour of modes entering the 
horizon at a time when radiation contributions cannot be 
neglected. During radiation domination (again assuming 
a small stress-energy of the vector field at these early 
times) and before recombination the coupled photon- 
baryon fluid will undergo oscillations of constant co- 
moving amplitude. On an expanding background, the 
associated behaviour of the metric perturbation will be 
to undergo damped oscillations, these in turn sourcing 
the vector field equation. Once again we find that the 
vector field tends to the growing mode solution. As for 
the longer wavelength modes, the initial conditions dic- 
tate that the scale dependence of Vb(fc) (and thus the 
shape of the power spectrum if V produces a dominant 
stress-energy perturbation) will be informed by the be- 
haviour of the metric perturbations. From (|32[) it can be 
seen that there will then be damped oscillatory terms in 
the vector field equation, there then, as in (|100[) . will be 
an integrated contribution to V with a functional depen- 
dence on the numbers Cj and n. Indeed we find a range 
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of parameters for which the damping of the metric po- 
tentials can become imprinted on Vo(k). If k 2 Vo(k) is a 
decreasing function of k, then the influence of the vector 
field will diminish on small scales. 

In summary then, modes of the baryon density per- 
turbation which enter the horizon while the metric per- 
turbations are being damped may diminish in size with 
increasing k at late times due to the effect of the damping 
of the metric perturbations on the evolution of the vector 
field, which, if yielding a realistic power spectrum, will at 
later times dominate the Poisson equation. Meanwhile, 
modes which enter the horizon after radiation may be ne- 
glected will grow as k? with increasing k. The demarka- 
tion between these regimes is the time of radiation-matter 
equality with an associated scale of the order k eq ~ r\~^ 
and we expect this to be observable in the resultant power 
spectrum. This is observable in figure ©) which plots 
numerical solutions for the matter power spectrum. 



B. Case 2 




k (h Mpc" 1 ) 



We now turn to the case where the initial conditions 
for E and V/rj are considerably greater than metric per- 
turbations. Recall that in the opposite case, we saw that 
the behaviour of the metric was instrumental in setting 
up the scale dependence of the growing vector mode. In 
this case, however, we will find that the eventual scale 
dependence of Vo(k) is predominantly determined by our 
initial values of E and V/rj. Again we first consider modes 
which lie outside the horizon during radiation domina- 
tion. The solution to V will again be described by equa- 
tion Again the decaying mode is neglected. Now 
though, we additionally may neglect the influence of the 
metric source. This imposes E = V/r) and: 

V = V( Vo ,k) (^-Y + (103) 

Clearly then, if we choose an initial condition 
A i (r]o, k) oc k~ q a dominant vector field will produce the 
following behaviour in 5(k) at late times for modes en- 
tering the horizon after radiation domination: 

s 2 (k, n) = f2(n)k 4 - 2 ^ (k« k eq ) (104) 

where the /^(f?) is determined by the solution to the 
field equations. 

Crucially, the solution (| 1 03[) in this case now even ap- 
proximately to smaller wavelength modes entering the 
horizon during radiation domination: the metric source 
terms remain but are not significant in altering Vo(k) 
and having suffered damping cannot appreciably influ- 
ence the solution at least until the vector field begins to 
act a a dominant source of stress energy. Therefore, the 
solution (|104|) is a good approximation across all modes 
of interest! In order to obtain a fit to the power spectrum 



FIG. 3: Density contrast squared as a function of wavenumber 
at z = 0. The dot-dash line is the ACDM prediction while 
the dotted line is the prediction for a baryonic universe with 
fib = 0.27. The solid and dashed lines represent baryonic 
universes f2(, = 0.04 with a dominant vector field for cases 1 
and 2 respectively. Data is from the Sloan Digital Sky Survey 



on the most observationally constrained scales (0.01-0.1 
Mpc -1 ) we must take q = —3/2 i.e. to have the same 
scale dependence as $ pr i m , and, by the fixed- norm con- 
straint, A . Immediately then we can see that in this 
case the resulting power spectrum will diverge on large 
scales (see figure (J3]))- 



VIII. COMPARISON TO TEVES 

There exists another relativistic theory which reduces 
to MOND in the quasistatic, weak-field limit: the TeVeS 
theory proposed by Jacob Bekenstein ([HI). As men- 
tioned, it has been shown [l5| that this theory can be re- 
written as a tensor- vector theory much like the model dis- 
cussed in this paper, the significant difference being that 
the coefficients Ci are now functions of A^A^ = A 2 - the 
non- vanishing vacuum expectation value of A 2 appearing 
not via a constraint but dynamically via the Ci(A 2 ) func- 
tions diverging as A 2 — > 0. It has been found that it is 
also the growth of the vector field perturbation V which 
is crucial in TeVeS for the formation of cosmic structure. 

Denoting the cosmological background value of A 2 as 
A 2 the field equation for V in TeVeS in a regime where 
a oc r\ m is: 
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V 



b -± V > + = S[<Z>,V] 

T) TJ 



where now: 



b 2 



m(3 - A 4 ) 

m(m — 1)(2 — A 4 ) 

+A 4 K B ~ 1 3m 2 (l - A 4 ) 



(105) 



(106) 
(107) 



where Kb is a constant of the theory. 

It was found in [191 ] that during the matter era the 
quantity A 2 is a function of time that will deviate very 
slightly from —1 towards as cosmic time increases. 
Again the metric source terms are initially subdominant 
during the matter era. If Kb is comparitively large then 
its contribution to (|105|l is suppressed. In this limit 
b\ — > 4 and 62 — * 2 leading to decaying homogeneous 
solutions for V: V oc r/ _1 and V oc r\~ 2 and the evolu- 
tion eventually being dominated by the particular solu- 
tion -*?7/3. 

Again we may choose the set of evolution equations 
where the contributions of the vector field to the evo- 
lution of other fields are characterized entirely by the 
energy density — a 2 5T° and anisotropic stress fc 2 (\& — $) 
(see [IDS) . These are given by: 



a 2 Sf° 



fc 2 (^-$) 



{I -A 4 ) 



A+ 



Hk 2 V 



K b ,2 



(1 - ,2 



A+ 



k 2 {2HV - 



2 
E) 



k 2 E (108) 
(109) 



Therefore so the system of evolution equations will be 
identical to (|79"]) with contributions from the vector field 
instead given by (|108[) . Recall that in the absence of 
the vector field, the baryonic contribution to the Poisson 
equation in the matter era is of the order S/r] 2 (where we 
have used that fact that 7i = 2/?/ during the matter era). 
The overdensity grows as n 2 and so the contribution is 
constant in time. In the event that the vector follows 
the particular solution proportional to r\ we see that the 
contribution to the Poisson equation proportional to E 
is also constant in time whereas the term proportional 
to HV grows very gradually due to the background be- 
haviour of A 2 . However, it was found that for modes 
entering the horizon prior to the end of radiation dom- 
ination that the value of VP appearing in the particular 
solution was sufficiently damped compared to the primor- 
dial value that the influence of the vector field remained 
too small to produce realistic growth of structure. There- 
fore, a more rapidly growing V is required. 

If Kb is indeed though, the contribution to 62 cannot 
be ignored. It acts to create a growing mode solution for 
V, with the exponent of 77 itself a function of a, eventually 
becoming greater than 1 and dominating the particular 
solution. It was found that this more rapid growth does 



allow for the growth of structure, indeed providing en- 
hanced growth over the ACDM model. 

However, it can be seen from the vector equation that 
the growing mode in the homogeneous solution only pos- 
sible if A 2 ^ — 1 . Therefore, the fixed-norm constraint in 
and of itself is not necessary for structure formation in 
a tensor-vector theory. A second noteable difference oc- 
curs in the rate of growth of the stress-energy of the vec- 
tor field. Recall that in the fixed-norm model, the time 
dependence of the background function Fk may signif- 
icantly alter the time evolution of the perturbed stress- 
energy whereas in TeVeS there is an analogous role played 
by (1 — A 4 ) I A 4 and a fixed number Kb , which itself must 
be small to allow a growing mode. 



IX. * - $ 

We can see from equation ([56")) that the difference in 
conformal Newtonian gauge potentials ('J — $) may be 
sourced by the transverse,traceless component of the vec- 
tor field's stress energy tensor. This component is pro- 
portional to the combination c\ + C3. Therefore we may 
expect the conformal Newtonian gauge potentials to only 
be equal if C\ + C3 = 0. 

Meanwhile, restricting ourselves to the parameter 
space where V tends to the growing monomial solution 
requires that the number 63 in equation (|43[) should be 
approximately zero. 

It can be seen from ([4"6")l though that 63 = is generi- 
cally satisfied if c\ + C3 = and C2 = 0. However, with 
this combination the vector field has no influence on the 
background evolution and thus no late time acceleration 
can occur. If we allow C2 to be nonzero, we can have 
63 = and ci + C3 = only when the number n is equal 
to 1/2. However, we have found that growing modes do 
not exist for this value in the matter era. Therefore, any 
growing monomial solution for V comes from a set of Ci 
such that C1+C3 7^ 0. Thus, from ([56"]) there will generally 
be an accompanying source for \P — $. 

This will lead the metric potentials $ and ^ to dif- 
fer in a scale dependent manner. This difference may be 
observable in physics at late times and at recombination. 
As discussed in ( [nj , [H , [13 ) at late times one may com- 
pare the deflection of light by large scale structure (an 
effect proportional to $ + \&) to the motions of galaxies 
in clusters (an effect proportional to $), thus deducing 

— <f> on large scales. Additionally, under the assump- 
tion of adiabaticity we have that the anisotropy in the 
CMB, AT(n)/T in a given direction n is given by 



AT(n) 
T 



1 



*(r?*,d*n)- / drj'(y + $)[?/, (770 - r}')n] 
J-n* 

( 11Q ) 

where 77* is the conformal time of last scattering, 770 is 
the conformal time today, and is the comoving radius 
of the surface of last scattering. 
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We can obtain a handle on the behaviour of the tem- 
perature anisotropy piO[) on scales which remain outside 
the horizon today. We first consider the term propor- 
tional to the value of the potential \& at recombination. 
Up to recombination, we assume that non-negligible con- 
tribution to the background energy density is due to the 
baryons and photons. In the absence of dark matter, 
radiation domination ends only after recombination. In 
analyzing the evolution of the potentials up to recombi- 
nation then we may neglect the baryonic contribution to 
the Einstein field equations. Therefore a oc r\. Further- 
more we will assume that we're at a time when the vector 
field has tended to the monomial solution for modes of 
interest and in doing so has not appreciably affected the 
other fields sourcing it (as one would expect from the 
early smallness of Tk) and so the contribution to the 
vector field in the system of field equations will be given 
by equations ([75]) and ([74")) (with m now equal to 1). We 
find this to be an accurate description of the evolution of 
superhorizon vector modes for the eventual production of 
realistic power spectra. As before, when considering su- 
perhorizon modes the velocities become decoupled from 
the field equations and terms explicitely proportional to 
positive powers of k are discarded, yielding: 



S' y = 4$' 



(111) 



3 _ 3 _ 

?7 jyZ 



3 
2?7 



a 5, + 3l s mv 3+p - en (112) 



where, again, the anisotropic stress due to the vector may 
be significant for these modes. Differentiating the second 
equation we may use the first to replace the derivative of 
the density contrast with a derivative of the potential, 
thus obtaining a second-order equation for <f>: 



$" + _$' = (5 + p - 6n)l s C(k)T 1 3+p ' 6n (113) 
V 

with solution 



=*<>(*) 



Is 



8 + p — 6n 



5+p— 6n 



(114) 



where $0 is the primordial amplitude of $ (it is identical 
to the integration constant as the vector field contribu- 
tion vanishes as rj — > 0) . 

It immediately follows (using equation ([74"]) ) that: 

V*) = Mk) + Is ( tt- 1 fi - l) ak)vl +p - 6n 
\8 + p — bn J 

Therefore the degree to which super-horizon modes de- 
viate from the primordial k dependence of the potentials 
at recombination depends upon the function For 
Case 1 initial data, it can be seen from ()101[) and (|76|) 
that for modes outside the horizon during radiation dom- 
ination, will have the same scale dependence as the 



primordial value of <E>. Assuming adiabadicity then we 
have £(fc) oc & pr im{k) oc fc~^. Thus in Case 1 the ef- 
fect of anisotropic stress on superhorizon modes will be a 
time dependent, scale invariant rescaling of the primor- 
dial value. As can be seen in (|110[) . a significant time 
dependence may produce a large scale CMB anisotropy 
today which is not approximately determined by primor- 
dial values (in contrast with the dark matter case). The 
effect is even more pronounced with Case 2 initial data 
where, by Q103P and ([75]) . we see that the k-dependence of 
£(fc) is approximately given by the initial k-dependence 
on E and V. It was found that a reasonable fit was ob- 
tainable for the linear power spectrum for V ~ fc~ 5 / 2 . 
However by (|1 14|) and (I115j) we see that this will lead to 
a correction to fc 3 <i> 2 ,fc 3 4> 2 , on superhorizon scales which 
diverges for small wavenumbers. 

We now consider the integral contribution to (|110|) . the 
so-called Integrated Sachs- Wolfe (ISW) [25[ . It may also 
be shown that, allowing for the different form of dom- 
inant background energy density, similar expressions to 
(|114[) and (| 1 1 5[) hold for superhorizon modes during the 
matter era, the sum of superhorizon potentials retaining 
a time dependence thus leading to a contribution to the 
integral. Recall that in ACDM cosmology, this integral 
becomes nonzero only after the onset of late time accel- 
eration. Again the Case 1 field emulates the initial scale 
dependence of the metric while for the Case 2 field the ef- 
fect of the anisotropic stress due to the vector field leads 
to a divergence on the largest scales. 

The latter behaviour is illustrated in figure ([U where 
potentials at z = are compared between the ACDM 
model and the Case 2 vector model depicted in figure 
([3"|). The wavenumber dependence of the potentials in 
the dark matter universe and Case 2 vector model agree 
for the range where the resulting baryon power spectrum 
is roughly similar. They are not of the same magnitude 
as in the Case 2 model the vector field has become the 
dominant energy density perturbation and so the matter 
perturbation is no longer simply related to the potential 
via the Poisson equation. For larger scales both poten- 
tials in the Case 2 model diverge from the dark matter 
case and from one another. As far as these initial condi- 
tions reflect the absence of turnover in the power spec- 
trum, they also lead to 'J, $, and their sum diverging on 
large scales, these effects growing with conformal time 
due to the induced time dependence of the potentials. 



X. DISCUSSION AND CONCLUSIONS 

We have shown that the vector field in the General- 
ized Einstein Aether theory can play the same role as 
in TeVeS and the same role as dark matter in sustaining 
the growth of structure during recombination. A growing 
mode for the vector field exists in the matter era for a 
wide variety of parameters and, through its associated en- 
ergy density and anisotropic stress, can become a domi- 
nant source term in the Poisson equation, thus sustaining 



13 



l 

0.1 k 
0.01 i 
0.001 i 
^ 0.0001 t 

% io-» 

10" 6 

io- 7 
io- e 

10" 9 
10 -io 

10" 11 



CO 




lO" 5 0.0001 0.001 0.01 
k (h Mpc" 1 ) 



FIG. 4: Metric potentials as a function of wavenumber at 
2 = 0. Solid and dotted lines represent fc 3 |$| 2 and fc 3 |*| 2 
each for a universe with (Qdm ~ 0.23,^6 = 0.04) (solid) and 
(£Idm = 0.0, fij, = 0.27) (dotted). Dashed and dot-dash lines 
represent fc 3 |>I'| 2 and fc 3 ^! 2 respectively for a universe with a 
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gravitational potentials which further act as sources for 
the evolution of the baryonic density contrast. The shape 
of the power spectrum is dictated predominantly by the 
k dependence of the growing mode which may either be 
determined by dominant initial data for the vector field 
(Case 2) or the manner in which an initially very small 
vector field reaches the growing mode solution via cou- 
pling to evolving gravitational potentials (Case 1). As in 
the case of dark matter, a turnover in the power spec- 
trum exists in the latter case due to the comparitive lack 
of damping of metric modes entering the horizon after 
radiation domination. 

Though the vector's stress energy tensor does contain 
terms proportional to the gravitational potentials (lead- 
ing to an effective time dependent rescaling of Newton's 
constant in Poisson's equation) we find that this is a sub- 
dominant effect. At late times and on large scales Pois- 
son's equation ([42]) takes the approximate form: 



V 2 $ 



4*Gp B + Y^K^n^V (116) 



where ps is the baryonic energy density. 

The dominant component of the right hand side at late 
times is found to be the term associated with the field 
V. In so far as the evolution of V is sourced by terms 
in the metric but not entirely fixed by it via constraints, 
its contribution to the Poisson equation is more akin to 
dark matter. 



It is interesting to contrast this with the weak-field 
quasistatic limit of the theory discussed in [l7|, wherein 
A M 's stress energy tensor is determined entirely by the 
metric and thus all observed mass discrepancies must be 
attributed to a modification of Poisson's equation which 
takes the form: 



V.(a*(|V$|/M)V*) = AttGpb 



(117) 



Therefore in this regime mass discrepancies can be in- 
terpreted as arising from a change in the relationship be- 
tween the dominant component of the matter energy den- 
sity (taken to be baryonic matter) and the gravitational 
field. In producing a realistic power spectrum, we find 
that the effective energy density of the growing mode V is 
considerable at the scale where overdensities become non- 
linear (k ~ 2hMpc~ l ). A similar situation exists in the 
TeVeS theory. This raises the issue, as yet unaddressed, 
of the degree to which V, instrumental in the growth of 
large scale structure, contributes to the mass of bound 
structure and, if so, whether it remains akin to a dark 
source in the Poisson equation (as in 111 6[) rather than a 
modification of how the gravitational field is sourced by 

In other words, covariant rc- 



baryonic matter (as in !117p 
alizations of the MOND program may well have uninten- 
tionally reintroduced dark matter, albeit non-particulate 
dark matter, via the back door. However, it is by no 
means obvious that the collective picture given by a uni- 
verse where a vector field gives rise to mass discrepancies 
would be degenerate with a particle dark matter model 
i.e the theory may resemble different types of dark mat- 
ter in different regimes but not a single type across all 
regimes. Indeed we have seen that the cosmological case 
involves significant anisotropic stresses at late times while 
the quasistatic case considered in [TtJ does not. 

We were able to clarify the analogy to the cosmological 
perturbations of Bekenstein's TeVeS theory. While the 
vector field is also nonvanishing in the cosmological back- 
ground in that case, it is seen that the growth of struc- 
ture requires the vector field be of non-fixed norm. As in 
the model discussed here, a growing mode in the vector 
field eventually dominates the evolution of the baryonic 
matter, its scale dependence being determined by the 
influence of metric source terms in its equation of mo- 
tion. However, in TeVeS the growing mode and vector 
field perturbed stress energy tensor will tend to have dif- 
ferent time dependences and so make for a potentially 
differing evolution of perturbations over cosmic time. In 
both theories the quantity ST° will generally retain a 
time dependence during the matter era, in contrast to 
the perturbed dark matter density in ACDM models. 

We now consider the parameter space explored in this 
paper. We have found that a power spectrum with a 
turnover and giving a reasonable fit to data tend to follow 
when the vector field and its time derivative are initially 
small. Following this, the degree to which the power spec- 
trum today can agrees with the data is dictated the val- 
ues of the parameters n and 7 (recall that T — j(—JC) n ) 
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and the Cj (see section fVII A[) . In the background we have 
seen (see figured]) through that some values of n and 7 in 
the modified Friedmann equation permit late time accel- 
eration. To first order in perturbations we have seen that 
the influence of n on the evolution of the matter and met- 
ric fields is via its combined influence on the exponent of 
the vector field growing mode and, along with 7, the time 
dependence of the quantity Fk to zeroth and first order 
which themselves act as sources in the Poisson equation 
and difference between Newtonian gauge potentials. In 
particular we were able to find values of n consistent with 
late time acceleration and a growing mode. We were un- 
able to produce a realistic power spectrum in the absence 
of dark matter for the value n = 1 (as hinted at by the 
lack of growing modes in the homogeneous solution to 
the vector equation for this value). 

The c,-, which determine the vector field's kinetic term, 
affect the degree to which the Friedmann equation is 
modified in the cosmological background. For instance, if 
ci + 3c 2 + c 3 = then there is no effect on the background 
evolution. To first order in perturbations, the Ci deter- 
mine the nature of growing mode. We have restricted 
ourselves to values of 63 (see equation (|46p) which led 
to an approximate monomial growth of the vector field. 
Deviations from this behaviour will have an integrated 
impact on the evolution of the other fields again, for in- 
stance, acting as a dark source in the Poisson equation 
but it is to be expected that different n, Cj will match 
observations of the matter power spectrum. 

It is worth noting that in this case as well as TeVeS a 
power-law growing mode emerges from a non-Maxwellian 
kinetic term (i.e. K. —F^F^a + c 3 7^ 0); it has 
been shown in [IH, [2y, [27| that such kinetic terms can 
arise from a change of variables at the level of the ac- 
tion. For instance in TeVeS the vector field kinetic term 
is Maxwellian when written as a bimetric theory but not 
when written as a single metric theory. Our choices for 
F(JC) were dictated by a desire for simplicity. It is to 
be hoped that a deeper grounding of the ideas discussed 
herein would fix the expected form of T . As all con- 
tributions of the vector field to the matter and metric 
evolution equations at 1st order in perturbations are pro- 
portional to the background value of Tk% the time depen- 
dence of this quantity can be expected to have a signifi- 
cant impact. We have seen that this will tend to produce 
enhanced growth for viable forms of a monomial J-(JC). 

Though we have chosen the sign of the term a = 
C\ + 3c2 + C3 to be consistent with subluminal propaga- 
tion of gravitational waves in the limit far from matter, 
in generic backgrounds the speed of gravitational waves 
will depend on F(JC) and it is not obvious that the the- 
ory would pass tests such as those discussed in [28| . Re- 
cent work [29, 3(| on the consequences of the theory with 
F oc K, have chosen the opposite sign of a where gravity 
propagates superluminally with respect to the preferred 
frame, thus avoiding any constraints from the above tests 
[3TI ] ; such a choice would change the requirements for a 



growing vector mode for more complicated F(JC) as well 
as leading to the background cosmology and weak field 
limit being described by the same sign tC, rendering them 
no longer independent regimes. 

Finally, the smallness of the mass scale appearing in 
the vector action M ~ cH n « Mpi is purely phe- 
nomenological. Though the value lends itself to modifi- 
cations of gravity associated with low gravitational fields 
in the weak field limit and at late times in the cosmolog- 
ical background (i.e. as H — ► M), it is fixed as a number 
in the action entirely by the former and lacks a deeper 
theoretical justification. A model wherein M itself is a 
function of cosmic time (through, say, a background time 
dependence on A 2 ) may have considerable effects on the 
background cosmology as well as time evolution of the 
perturbed vector fields stress energy tensor. 

As the size of perturbations is constrained on large 
scales by the baryon power spectrum at late times and 
the CMB anisotropy at recombination (allowing for the 
ISW effect), a full simulataneous modelling of each in the 
presence of the vector field should allow further discrimi- 
nation between models and the degree to which they can 
give an account of mass discrepancies compatible with 
data in precision cosmology. 

In summary we have considered a set of models of a 
universe with a timelike vector field with a noncanonical 
kinetic term. We have found that the model on one hand 
exhibits behaviour typical of 'modified gravity' theories 
in high symmetry (for instance modifying the Friedmann 
equation in the FRW background and the Poisson equa- 
tion in the static spherically symmetric case [13] )• Addi- 
tionally we have found that the extra degrees of freedom 
introduced through covariant realization of the former 
behaviour behave in a manner more akin to dark matter 
in perturbations around the FRW background. In par- 
ticular, we find that at late times the stress energy tensor 
components of the field V can become the dominant com- 
ponent in the cosmological Poisson equation. It is not 
clear, for instance, that this dominant component should 
be found to align with displacement of the subdominant 
baryonic matter in large scale structure. Even in so far as 
the field V may resemble dark matter in the linear regime 
we have found a number of aspects in which its account 
differs from that of cold dark matter. We hope that this 
may help point the way towards a general approach to 
observationally distinguishing between the concordance 
model of cosmology and modified gravity. 
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